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SUMMARY 

Fluid forces on a centrifugal impeller, whose rotating axis whirls with 
a constant speed, have been calculated by using unsteady potential theory. 

To simplify the problem, it is assumed that the flow is incompressible and two 
dimensional, the impeller is surrounded by vaneless diffuser and the 
eccentricity is constant and small. Calculations are performed for various 
values of whirl speed, number of impeller blades and angle of blades. Some 
typical examples of results are shown in vector diagrams. Although the results 
of this paper are obtained for simplified geometries under idealized 
assumptions , they will suggest us the nature of fluid forces acting on whirling 
impeller. 


INTRODUCTION 

Recent advances in rotor dynamics made it possible to analyse very 
complicated rotor system like multi-stage pump rotor and multi-bearing 
turbogenerator shafting. In these analyses the largest weakness seems still to 
lie in the assumption of radial fluid forces acting on rotor elements such as 
seal rings, balance pistons and impellers. Many studies have been made and are 
being carried on this subject; just to mention a few, contribution of Black to 
pressure seal and centrifugal impeller (ref. 1,2) and of Acosta to centrifugal 
impeller (ref. 3). 

. As to radial forces acting on centrifugal impellers, numerous studies and 
informations have been accumulated on their steady part caused mainly by asym- 
confi S ura tion of volute casings. To the contrary limited studies (ref. 
2,3) have been published on the unsteady part of radial forces, which is caused 
by the whirling motion of rotating impeller. This is the very information 
sought for the improvement of rotor analysis. Most of them are, however, 
quasi-steady analyses, in which the effect of shed vortices is disregarded. 

In this paper, the radial force of a rotating impeller with arbitrary 
whirling motion is reduced from the flow fields in and outside of the impeller, 
applying familiar singularity method for the calculation of unsteady potential 
flow. Though the flow in a real centrifugal turbomachine is very complex 
because of asymmetric geometry of casing and strong viscous effect in impeller 
passage for example, following simplifications and assumptions are introduced 
in the present theory. 


*Not presented at workshop. 
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The fluid is incompressible and inviscid, and the flow is two 
dimentional and irrotational in the absolute coordinate system. 

The region outside of the rotor is unbounded and free of vanes. 

The ratio of whirl eccentricity to the radius of the impeller is so 
small that its square can be neglected. 

The flow rate, the prerotation of suction flow, the rotational speed 
of the impeller and the whirling speed of the shaft center are 
constant . 

Blades of impeller have no thickness and have the shape of 


logarithmic spiral. 

(6) Steady part of the flow satisfies the condition of shockless entry 
at the leading edges of the blades and Kutta's condition at the 
trailing edges is assumed. 

(7) Free vortices are shed from the trailing edges and are carried 
downstream along steady streamlines with steady velocities. 

After having derived formulations of radial forces for given parameters 
of whirl to rotational speed ratio, number of impeller vanes and blade angle, 
typical numerical results are demonstrated in vector diagrams. 


SYMBOLS 


N 


u) 


n 


e 

r, 

Y 

i 

j 


number of blades 8 angle between radius and blades 

inside radius of impeller r 2 outside radius of impeller 

angular speed of impeller, positive clockwise 
angular speed of whirling motion, positive clockwise 
eccentric radius Q 

prerotation P 

strength of vortices z 

complex unit with respect to coordinate 
complex unit with respect to time 


flow rate 
density of fluid 
complex conjugate of z 


Subscripts 

n component normal to blade t component tangential to blade 

s steady component 0 quasi- steady component 

w component caused by free vortices 

co component caused by rotation effect 
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ANALYSIS 


Basic Consideration 


wit-h ceatrifugal im Peller shown in figure 1 is considered. It rotates 

radius e°and a a conft r SP6e ? “ ^ US C6nter Whirls with a constant eccentric 

velocity field U ^ ■ angUlar Speed The bounda ^ condition to the 

shouldL^IJr a l hlS case unstead y^ and the flow around the impeller 
nould be determined by unsteady analysis. 

Let an arbitrary point on a reference blade be A. Then the geometry of 
the points A, 0 and O' are shown in figure 2, where point 0 and O' are the 
center of rmpe n er and whirling motion, respectively. From the viewpoint of 
flow field, we interpret that the inlet flow to the impeller is unaffected by 
the whirling motion of the impeller, that is, the location of source Q (repre- 
senting flow rate) and circulation T p (representing prerotation) is fixed on 
point 0 . Defining the angles 3, 0 1 and 02 as shown in figure 2, and 
introducing the lengths AO •=*, AO=r and 00' =e, £ is expressed as 

1 = t r2+ 2 r cos(0 2 - 0 1 )+ e 2 ] 1/2 = r [1 + e/r cos(9 2 - 0 )] 

The velocity at A induced by source Q in the direction of O'A is 


v = Q/ 

Q / 2 TT l 


7 2 71 r [1 ' % COS(0 ' 


6 ,)] 


Similarly, the velocity vp at A induced by circulation T is 
r P 


v r '■ Pp/ 2 „ r t> - V r cos(e 2 - 


ep] 


0 ) 


( 2 ) 


where its direction is normal to the line O'A (positive clockwise). 

The velocities vq and v p are divided into two components respectively, 
ihey are y Qn and v pn which are normal to blade, and v 0t and v which are 
tangential to blade. As the reference blade intersects the line OA at angle 6 
e intersection angle between blade and line O'A is given approximately by 
R + rrAO -0 1 Using the above re i a tions, v Qn> v^, v^ and v^ can be 


3 + (e/r) sin(0 -0 ) . 
expressed as 1 


Qn 


•= Q 


2 tt r f sin B + % sin(0 2 - g)] 


V rn ' P / 2 ^ r fcos B - / r cos(0 2 - 0 r g)] 

2 tt r t cos 3 - £ / r cos(0 2 - 0 1 - g)] 
P/ 2 tt r f sin 3 + £ / r sin(0 2 - 0 1 _ g)] 


v = Q/ 

v Qt ‘ ' 


Tt • 


(3) 

(4) 

(5) 

( 6 ) 


When point A on the blade is observed from a fixed coordinate system A is 
moving in the normal direction to OA with velocity ru>, and also in the normal 
irection to 00 with velocity eft. Both velocities are considered positive 
c ockwise. On the other hand, when observed from the relative coordinate 
system fixed to the rotating impeller, absolutely still field is considered to 
have velocities -rw and -eft in the corresponding directions. These are divided 
into two components m the same way as vq and v r . Using the relations 
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d0 1 /dt=-n and de 2 /dt = -aj, normal and tangential components and v^ are 


v = -r a) cos g - e ft cos(0_- 0,+ g) ( 7 ) 

(jjn 2. I 


v = r a) sin g + e sin(0„- 0,+ g) (8) 

wt 2 1 

Hereafter only relative coordinate system fixed to the impeller will be 
used. Every vector quantity such as location vector 0A, force and velocity is 
indicated in this coordinate system by a complex number, whose real and 
imaginary (with imaginary unit i) part represents its component in two 

orthogonal axes. . 

Imaginary unit with respect to time, j, is also introduced to indicate the 

phase difference of oscillating quantities. Since two complex units l and j 
have completely different meanings, they should be distinguished carefully in 
the reduction of equations. 

Let the location of A be z in the complex plane, normal component vq^, 
and v wn of equations (3) , (5) and (7) can be written as 


V Qn ^2 tt |z| 
v I-* =r 


[sin g 


j £ z , j 3 j t. 

+J / , 2 e J e JV J - = v- ,+ 


e z 

P -*p/t „ , i [cos B - / .2 e 

fn V 2 tt \z 1 7 


j Qnl Qn2 
jB jCu-flJt, 

e ] j v fnl V fn2 


e z 


v = [- |z| co cos 3 - k i Q e J e 
con L 1 1 |z| 


-jBJ (w-fi)t 


1 . =v + v 0 
J j coni con2 


(9) 

( 10 ) 
( 11 ) 


Since only absolute values of above three quantities are of interest, they are 
expressed as a complex with respect only to j. Further, tangential component 
v Qt , v r and v wt of equations (4), (6) and (8) are written as 


'Qt * 'l *1*1 [C ° S 6 


£ Z / , pjSpl (W-fl)t, rZ, 

\z\ 2 e e J j 1 M 


iB, 

. e ] . 
z 1 


=V Qtl + V Qt2 


(12) 


r . , . „ i e z, iB j(w-f2)t 1 r z. iBi 

v rt = - p/ 2 TT ,z| [sin B + 7 |Z| 2 e e J ^’ [ / - e 


'j 1 |Z| i 


_v rti + v rt2 


(13) 


v wt = t l z l w sin 3 + e ft j /| e 


-jBJCu-fiJt, ,Z, iBn 
e J j 1 \z\ e J i 


= V * 1 + V . ~ 

tot 1 wt 2 


(14) 
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All of the equations from (9) to (14) consist of the steady term indicated by 
subscript 1, and the unsteady term by subscript 2. Unsteady components are 
first order quantities in comparison with the steady component. 

Since all the unsteady terms contain this problem can be thought 

to be the periodic phenomenon with angular speed (w-ft) . All quantities within 
[ are complex with respect to j, and therefore their argument corresponds 
to the phase of sinusoidally oscillating phenomenon. 

Steady Flow 

As thickness of blades is disregarded, blades can be replaced by the 
distribution of vortices. The distribution of steady vortices y s has to 
satisfy the boundary condition that there is no crossing flow through blades. 
This leads to 


,Z T M N 

^ u,) Re( 7T 


iB 


|z| 


3 dZ ' + V Qn 1 + v mi + V ojn 1 = 0 ( 15 ) 


where and Zj are the complex coordinates at leading and trailing edges of 
the blade. The integral should be carried out along the blade. For the 
calculation, of the integrand of the above integral, the following relation is 
also used 


£ 


z 


1 

Z 1 e i 2fTk/N 



N - 1 
z 



(16) 


Combining equation (15) and the Kutta's condition at the trailing edge, y s can 
be determined. As shockless entry at the leading edge is assumed in this 
analysis, Ys is restricted to finite value. Once y s is obtained, the velocity 
v(z) at an arbitrary point z can be calculated by 


v ( z) 



C 2 ' ) 


2 TT 


N- 1 
z 


( 


N 

z - 



d z ' + 


2 (Q- i Tp) 
2 tt |z| 2 


+ i 


z (J 


(17) 


Unsteady Flow 

First we consider the flow field shown in figure 3, in which N discrete 
vortices with phase difference 2 tt/N are ranged on a circumference with equal 
distance. Let the location of a reference vortex be z', the complex conjugate 
velocity at an arbitrary point z is given as 


f(z,z') = u -i v= 


2 TT 


£ 

k = 1 


- j 2 it k/ N 


z - z 


i2irk/N 
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(18) 


1NT 

4ttz‘- 


(z/z* ) 


N/2-1 


(z/z ' ) N// “ 


+ (z’/z) 
- (z'/z) 


N/2-1 
N/2 ' 


, / n ^ 2 - 1 

•(z/z') 


(z/z')^ 


(z*/z) 

(z’/z) 


N/2-1 

N/2~ 


] 


The above function f(z,z’) is very useful for determing unsteady velocity, 
because the unsteady distribution of vortices y u has also the phase difference 
2 tt/N at the corresponding points on the blades. 

Unsteady distribution of vortices y u consists of quasi-steady term y 0 and 
wake effect term y , 


Y u = Y 0 


y. 


The induced velocity by Yq ^as to cancel 
(10) and (11) on each blade, 
equation as 2 

f \(z.) Im[ f ( z , z ' ) e lp, z/ |z| ] 


Qn2 ’ 


(19) 

f vn2 anc * v u)n2 equation 

i ^ 1 1 _ „ 1 


This relation is expressed by a integral 
IB 


dz’ 


v Qn2 + 


V Tn2 + 


(jon2 


0 


( 9 ), 


(20) 


Introducing the Kutta's condition at the trailing edge, y Q can be finally 
determined . 

The free vortices are assumed to be carried downstream along the 
steady streamline starting from the trailing edge, which can be calculated by 
equation (17). All wake streamlines have the same geometry, but free vortices 
on them have also 2 tt/N phase difference with one another. Considering the fact 
that the strength of free vorticies is inversely proportional to the velocity, 
the distribution of free vortices y^ of the reference blade can be written as 


A = c 


v(z T ) 

v(z) 


e j (w-fl)t e -j h(s) 


( 21 ) 


Here, C is a complex constant with respect to j, |v(z)| the absolute velocity 
on wake streamline, s the length of the streamline measured from the trailing 
edge z T . The function h(s) in the above equation is defined further as 


h(s) 



ds 1 / | v ( z ’ ) | 


( 22 ) 


Assuming C=1 provisionally, the velocities on the blade induced by free 
vortices, v wn are written as 


v 

wn 


(z) 


= Im[ 



v(z T ) 
v(z' ) 


h(s) 


f(z. 


Z ' ) d Z ' 


i6 


Z / |Z| ] 


(23) 


To cancel this normal component on the blade, the unsteady distribution of 
vortices y wl (y w i=y w /C) is to satisfy the following integral equation 


f \,P z, ) In| [f(z,z') e lB z/ |z| ]dz' + v (z) 
j 7 w 1 wn 

J L 


0 


(24) 


I'he above equation should be solved under the additional condition that the 
vortex strength is continuous at the trailing edge. 
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Introducing the circulations T s , Tq and r ,, which are caused by y«> YO 
an d Ywl> the total circutation T is 


r 


r + 

s 


r + c r 


vcl 


(25) 


Differentiating both sides of the equation (25) with respect to time and using 
d/dt=j (cu-fi) , 


dr/dt = j O - fl) (r 0 * c r ul , (26) 

is obtained. Further, Kelvin’s law for the conservation of vortices leads to 
the relation that the rate of total circulation change is equal to the total 
shed vortices from the trailing edge within unit time. This is expressed by 

- dr/dt = C|v(z T )| e'j h(0) (2?) 

where |v(z T )| is the steady velocity at tlie trailing edge. Substitution of 
equation (27) into equation (26) yields, 

-j (w - m r 0 

C = 1 Cu> - si) r wl +| V (z T )| t 28 ) 

Once C is determined in the above, whole unsteady flow field can be solved 
completely. 


Forces Acting on the Impeller 

Forces acting on the impeller can be determined by integrating the 
pressure distribution on the blades. The assumption of infinitely thin blades 
makes this calculation comparatively simple. Let A and B be upper and lower 
surface of a blade point, the pressure difference between these points is 
expressed as 

p B /p - p A /p -(V- V )/2 * 3 »a- *„>/3t 


= ( V A + V R ) Y/2 + j(u)-fi) 




Cy o + V dz 


(29) 


where v A and v B are the local velocities, and <|> A and <|) B are the velocity 
potentials. Here (v A +v B )/2, the mean velocity of the both sides, is denoted 
y m and is separated further into steady component U ml and unsteady component 


U, 


'm2- ^ml can be obtained from equation (17), 


ml Z A 3 ' V < Z A > 


(30) 


1121^113^4 nlf; ed u y adding v Qt2( z A)> v rt2( z A) and v Wt2 (z A 
(12), (13) and (14) to the tangential component of the velocity 
Y 0 . Y w and Yj . 


from equation 
induced by 
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2 (z A } = [[// (Y 0 (z’)n w (z*))Re[f(z A ,z')e lB z A /|z A !].dz' 

/ °0 . - 

Y 1 (z')Re[f(z A ,z')e l8 z A /| 2 A |] i dz-] e j(w - Q)t . • [z A e iB / | z A | ] . 
Zrp - 3 


* v Qt2 lz A 5 * v rt2 (z A ) * Vi'Y 


Consequently, the first term of the right side of equation (29) may be written 


(V A , Vg) V/2 = (U ml * U m2 ) CY S * Y 0 s V 


U ,y„ + U , (Y n + Y ) + U „Y 
ml's ml 0 ' w' m2 s 


where U m 2 (Yn + Y )> which is second order quantities, is neglected. Integrating 
the pressure difference along all blades, and adding them up as vector 
quantities, the total force acting on the impeller is given as 



Obviously the total force obtained is purely oscillatory without any steady 
component. This is because steady component vanishes by adding up for all 
blades . 

It is well known from thin airfoil theory that the direction of lirt 
acting on an airfoil is perpendicular to the inlet flow direction, and the lift 
cannot be calculated by merely integrating the pressure difference between 
upper and lower surfaces along the airfoil. It is because the infinite 
strength of vortex at the leading edge generates a finite force in the 
chordwise direction. In this analysis, y s does not become infinite at the 
leading edge because shockless entry is assumed throughout. Even in this 
case/ y 0 and y w become infinite at the leading edge. Calculation shows, how- 
that the effect remains always within second order quantities, and thus 
can be neglected. From the above considerations, it is evident that the 
integration of the pressure difference leads to the force acting on the 
impeller. 
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Method of Numerical Calculation 


In the practical calculation the vortex distrbution on blades and wakes 
are replaced by finite number of discrete vortices, thus transforming the 
integral equations to simultaneous equations. In the calculation of unsteady 
vortices, it is noted that the coefficients of simultaneous equations become 

r6SpeCt Numerical integrations should resort to trapezoidal 

method by separating the integrands into real and imaginary part The intP 

fro^T °f- Yl 6quati0n and ( 31 ^ can be terminated at a proper distance 

rom trailing edge, since the effect of free vortices diminishes rapidly as 
they get away from the impeller. y y 


Numerical Results 

Calculations are carried out for blade angle 6 of 30°, 45° and 60° 

Some examples are given in figure 4,5 and 6. Figure 4 is the result for the 
case of N=6, r 1 /r 2 =0.5, $=45°, r p =0 and the flow coefficient ip=0. 25. The 
absolute values and the directions of dimensionless forces acting on the 
impeller at the moment when 00' is in the direction of x axis are shown for 
various values of fl/w. Though the forces acting on each blade are unsteady, 
the combined forces become always constant, provided the number of blades is 
greater than three. Observing from the relative coordinate system, the force 
changes its direction clockwise with the angular speed (ft-w) . While, observing 
from the fixed coordinate system, the direction changes clockwise with the 
angular speed Q. The angle between force vector and 00' remains constant. 

=30°, r p =0 and if>=0.443, 
=0 and $=0.1443. 

. . ynamic damping force 

acting on each element of the rotor plays a predominant role. In the present 
case damping force corresponds to the component of fluid force on the impeller, 
which is parallel and is in opposite direction to the whirling velocity of 
shaft center. When damping force is negative, that is, when the fluid force 
has such component that pushes the shaft center toward its whirling motion, 
there exists a danger that the rotor system runs into severe whirling. 

In tlie results shown in figures 4, 5 and 6, damping force corresponds to 
the fluid forces in +y direction for fi>0 and in -y direction for £K0. As seen 
from figures, fluid force acts as damping throughout the calculated cases. 

When Q, is equal to u)(f2/io=l), flow through the impeller becomes steady and 
no free vortices shed out from the blades. In all other cases (ft/co^l) , blades 
are always under unsteady condition and the influence of 3<f>/3t term in equation 
(33) cannot be neglected. Quasi-steady treatment seems to be inappropriate 
especially for centrifugal impeller, because the wake vortices shed from one 
blade pass in the proximity of the next blade, thus resulting in a substantial 
interference. This tendency becomes more remarkable when blade angle is large. 

Besides the examples shown in figures, calculations were performed for 
various combinations of parameters. From these results it becomes clear that 
the number of blades has little influence on the total forces. 

CONCLUDING REMARKS 

Fluid forces on a centrifugal impeller, whose rotating axis whirls with a 
constant speed, was solved by the singularity method. Some typical examples of 


Mgure b is the result for the case of N=6, r 1 /r 2 =0.5, £ 
while figure 6 for the case of N=6, r 1 /r 2 =0.5, 6=60°,’ IV 
In the stability analysis of a rotor system, hvdroc 
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results were shown in vector diagrams. Although the results of this paper are 
obtained for simplified geometries under idealized assumptions, they will sug 
pest us the nature of fluid forces acting on whirling impeller. 

The present study is expected to be extended to those cases in which the 
condition of shockless entry is eliminated and the effect of volute or guide 
vanes is taken into account. Experimental study is also under way by the au- 
thors and is expected to furnish data to be compared with the present analysis. 

This study is supported by Scientific Research Fund of the Ministry of 
Education and by Hitachi, Ltd. 
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Fig. 1 Schematic geometry of a rotating 
impeller with whirling motion 
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Fig. 2 Relation between center of 
impeller 0, center of whirling O' 
and blade point A 


Fig. 3 Discrete vortices with 
phase difference 2tt/N ranged on a 
circumference with equal distance 



Fig. 4 Vector diagram of radial forces on whirling 
impeller (N=6, 6=45°, r 1 /r ? =0.5, r p =0, ^=0.25) 
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g. 5 Vector diagram of radial forces on whirling 
impeller (N=6, 6=30°, r^r^O.5, r p =0, ip=0.433) 


F 



g. 6 Vector diagram of radial forces on whirling 
impeller (N=6, B=60°, r./r =0.5, r p =0, ip=0.1443) 



